Antilocalization in a 2D Electron Gas in a Random Magnetic Field 
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We construct a supersymmetric field theory for the problem of a two-dimensional electron gas in 
a random, static magnetic field. We find a new term in the free energy, additional to those present 
in the conventional unitary cr-model, whose presence relies on the long-range nature of the disorder 
correlations. Under a perturbative renormalization group analysis of the free energy, the new term 
contributes to the scaling function at one-loop order and leads to antilocalization. 



PACS numbers; 72.15.Rn,73.20Fz,73.23.-b 

A question of long-standing controversy has been the 
influence of a random, static magnetic field on a two- 
dimensional electron gas. This problem is relevant to a 
wide range of physical situations, including the quan- 
tum Hall effect in the half-filled Landau level ^ and 
high-Tc superconductivity in the gauge-field description 

. Following the application of a succession of analyti- 
cal methods to this problem (e.g. |^), Aronov et al. ^] 
constructed a cr-model for 2D electron sea in a delta- 
correlated random magnetic field (RMF). They found 
that the a-model was identical to the one derived previ- 
ously for disordered systems with broken time-reversal in- 
variance (unitary ensemble) and therefore concluded 
that all states are localized as for the unitary ensemble. 
In particular, they did not find the additional term in the 
CT- model proposed in Ref . ; this term appeared in Ref. 

as a result of an incorrect relation between a corre- 
lator of Hall conductivities and the average longitudinal 
conductivity. Although recently some arguments in fa- 
vor of delocalization by the RMF were presented , it 
is apparently commonly believed that all states of a 2D 
electron gas in the RMF are localized and Ref. gives 
a proper description of the model. 

In this paper, we reconsider the problem presenting a 
more careful derivation of a proper cr-model. It repre- 
sents a generalization of the supermatrix cr-model |^ to 
the case of a disorder potential with long-range spatial 
correlations. For the RMF model, a new term (differ- 
ent from the one of Ref. |^) appears in the free energy, 
additional to those present in the conventional unitary cr- 
model. The presence of this term is a consequence of the 
breaking of a "gauge" symmetry due to an anomaly at 
large momenta. The long-range nature of the vector po- 
tential correlations corresponding to the RMF are crucial 
for our derivation. Under a perturbative renormalization 
group (RG) analysis of the free energy, the new term 
contributes to the scaling function at one-loop order and 
results in the antilocalizing behaviour. 

A related model which displays electron delocalization 
is the random flux model which describes, for example, 
electron hopping on a bipartite lattice structure with link 



disorder (see e.g. ||,^). Here a tendency towards delo- 
calization is displayed as the band center is approached 
[|l^,0, due to the existence of a chiral symmetry at the 
band center. However investigations of a lattice formula- 
tion of the RMF problem (equivalent to the random flux 
model at energies far from the band center) remain con- 
troversial: conclusions are divided between the localiza- 
tion of all states [0 and the existence of a critical region 
|ll| . The existence of a metal- insulator transition for 
an RMF model with a spatially-correlated magnetic field 
has also been asserted recently . It is also of interest 
that, even for the case of one dimension, recent analytical 
|l3| and numerical ||lj] work supports the existence of a 
metal-insulator transition in the presence of sufficiently 
long-ranged disorder correlations. It is well-known that 
spin-orbit impurities in 2D also lead to antilocalization 
|l5| , although by an entirely different mechanism than 
that considered here. 

In the following we consider the Hamiltonian, 

n= {-iV - eA/c)^ /{2m) ~ ep. (1) 

A delta-correlated magnetic field leads to a vector po- 
tential, A, with a long-ranged and transversal correlator, 
{A\t)A'{v')) = 2{mc/efV'^{v - r') where 
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for some strength 7. Here ep and vp are the Fermi- 
energy and Fermi-velocity and k <C 1 is a cutoff that 
renders finite the otherwise-infinite range of the disorder 
correlations. At the end of calculations one should take 
the limit k — > 0. This cutoff appears |3|j^ in the di- 
vergent single-particle lifetime: for example, the simple 
Born approximation yields = 2€.p^ / n while the self- 
consistent Born approximation (SCBA) yields the less di- 
vergent T^^g^ — 2tp{'^\n{\ / k) / ttY/'^ . The inverse trans- 
port time however remains convergent, = ep"f. We 
emphasize that our method remains valid for an arbitrary 
choice of the correlator V^^ (r) , with minor modifications 
for scalar rather than vector potential disorder. 

We now derive the free energy functional. Following 
standard methods and notation H , an averaged product 
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of Green's functions may be expressed as a functional in- 
tegral weighted by a Lagrangian that is quadratic in an 
eight-component, supersymmetric i/i-field, 

C^i f ^(r) (no - — rsA.V - -^aA ^{v)d'v. (3) 



Here Ho = e+eF + V^/ {2m)~ujA/2, and and A are the 
z-Pauli matrix in time-reversal and retarded-advanced 
spaces. We average over A (neglecting the term in as 
it gives only a slowly varying contribution to the chem- 
ical potential) to induce a term in the Lagrangian that 
is quartic in the ip-fields and, in contrast to the case of 
short-range impurities, non-local in position. Decoupling 
this term via an 8x8 matrix Q(r, r') and integrating over 
the ■(/'-fields, we come to the Lagrangian 
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Str\n(iHo + Vr,r'Q{r,r'] 



+ istr(Q(r,r')V;,r'Q(r',r) 
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(4) 



with K.r' = -(l/2)^^^.l/»^(r-r')(V;-V;,)(V^r-Vi:,). 
Also Q = Q\\ + iQ±, where Q = Q|| + Q± and Q|| 
(Q-l) commutes (anticommutes) with t^. The Q-matrix 
satisfies the standard symmetries Q{r,r') — Q{r',r) — 
KQ+{y',y)K, with K defined as in Ref. §. 

To search for the saddle-point to Eq. (|4|), we take 
Q± = and Q to depend on r — r' only. After Fourier 
transforming, the saddle-point equation reads Qp = ^p, 
where 
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and ^p = p^/(2m) — ep- Writing Qp — Ap, the required 
form for Ap is in general complicated. The analysis sim- 
plifies however in the limiting cases of k 3> {epT)^^ and 
K ^ (e/T-r)"^, where r is the mean free time that is to 
be found from the solution. The latter case is relevant 
here. In both cases we may take Ap = i/ {£_p + iA/ (2t)), 
while keeping epT ^ 1. In the limit k <C {epT)~^ one 
comes to the SCBA solution for the scattering lifetime, 

T — TSCBA- 

To consider fluctuations around the saddle-point, we 
separate them into three types: (a) hard massive, associ- 
ated with the fluctuations of the eigenvalues of the Q and 
the cooperons, with the mass that is very large in the 
limit K — > 0, (b) soft massive, with the mass r^7^ that re- 
mains finite, and (c) massless. The condition epT ^ 1 
ensures that fluctuations of type (a) are irrelevant and so 
Q^{r,r') = 1 and [Q,r3] = 0. To obtain a functional for 
types (b) and (c), we write 



Qp(R) = U (R) g(°)c7 (R) ; Q(°) = K.(R)ApK (R) 



with R = (r + r')/2 and n = p/\p\. In Eq. (|), U (R) 
is independent of n, and Qp ^ (R) contains only non-zero 
harmonics in n around the Fermi surface. The superma- 
trix q'^^ (R) can be parametrized, for example, as 

Q(,") (R) - Ap (1 + zP„ (R)) (1 - zP„ (R))-i (6) 

with / dnPn = 0, P„ = -P„, {Pn,A} = 0. Both [/ (R) 
and Pn (R) vary slowly with R, that is, on length scales 
longer than / = vpT. 

The parametrization for Qp (R) guarantees that the fi- 
nal free energy functional describing fluctuations will be 
invariant under global rotations U (R) UoU (R) with 
Uq = Uq^ independent on coordinates, this invariance 
being present in the initial Lagrangian, Eq. (^). 

As we are interested in the low-frequency behavior of 
transport coefficients the massless modes are of greatest 
interest for us. However, we cannot simply neglect the 
soft massive modes (b), and our aim now is to integrate 
them out, thus reducing the free energy to an effective 
functional containing only U (R) . Due to the soft mass 
of the (b) modes it is sufficient to apply a Gaussian ap- 
proximation in Pn (R) . Higher order terms give a small 
contribution provided the inequality epTtr S> 1 is ful- 
filled. 

The crucial feature of this integration is that, as a con- 
sequence of the long-range correlations of the disorder, 
Eq. (^), it requires a cutoff at large momentum k (short 
distance) which remains in the final functional. There is 
some freedom in how to choose the cutoff for fluctuations 
of Pn (R) as the symmetry under the global rotations is 
already guaranteed by the choice of parametrization and 
there are no other symmetries to be respected in the ini- 
tial Lagrangians, Eq. (|^, The simplest way is to write 
this cutoff as the distance where Itr = vpTtr Iq 1-, 
below which the fluctuations are suppressed 



Pnk = for |k| > l^' 



(7) 
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where P„k are Fourier components of Pn (R) = 
X^k exp(ik.Rr3). As we have to calculate only Gaus- 
sian integrals over Pnk, Eq. (^ is very convenient for 
explicit computation. At the end of the calculation we 
will take lo —>■ I- 

We remark that to recover the conventional cr-model 
for short-range (delta-correlated) impurities we may 
set Vn = I in Eq. (||) so that Q is a function of R 
only. It then follows that the free energy possesses a 
local "gauge" invariance (LGI) under the transformation 
U{R) -> [/(R)/i(R), where [/i(R), A] (note the term 
"gauge" is used here in a separate sense from that of 
the original vector potential). It is the loss of this in- 
variance in the case of long-ranged disorder that leads 
finally to antilocalization. Before explaining this point 
we present the expression for the free energy functional 
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F[U] obtained after integration over the non-zero har- 
monics Pnk^ 



F[U] 



— Str [D{WQ (r))2 + 2iujAQ (r) 
8 
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(8) 



where D = u|,rtr/2 is the classical diffusion coefficient, 
Q(r)^U (r) AU (r), and (3 = nr/ (2r) > 1. 

The first line in Eq. corresponds to the conven- 
tional cr-model for the unitary ensemble that predicts 
localization. This part has been derived in Ref. M for 
the model in question. The second line in Eq. (^) is 
the main result of the present paper. The first line re- 
spects the LGI but the term in the second line violates 
this invariance and leads eventually to antilocalization. 

At first glance, the violation of the LGI seems im- 
possible because, on replacing C/(R) —> C/(R)/i(R) in 
Eqs. (|^, H), we can change the variables of integration 
Pn (R) h (R) Pn (R) h (R) when integrating over the 
non-zero harmonics. This would remove h (R) from the 
final expression for F[U]. However, this change of the 
variables would replace Eq. (j^) by a more complicated 
expression for the cutoff. The LGI is therefore guaran- 
teed only if the cutoff is not important, that is, when no 
ultraviolet divergencies appear in the integral over the 
non-zero harmonics. 

Now let us explain the most important steps of the 
derivation of Eq. (^. First, we perform an expansion 
of the logarithm in the Lagrangian (Q), in both TnVTn 
(where Tn = UVn) and the collision integral (devia- 
tion of T4,r'Q(i", I"') from its saddle-point value). If we 
keep to first order in both of these terms, we recover 
the Muzykantskii and Khmelnitskii a-model with ballis- 
tic disorder . Although such an expansion is sufficient 
to describe a long-ranged potential (for which rtr 3> t), 
for the general case including a short-range potential we 
require the corrections to the free energy in Ref. |l|] 
obtained by continuing this expansion to second-order in 
both terms. The resulting terms in the free energy then 
complete the list of those that contribute to Gaussian 
order in Pn. 

We find it convenient to Fourier transform from Pnk 
to angular harmonics, Pnk = X^m *5xp (iTi^ra). Ex- 
pressing the collision integral in terms of these harmonics, 
we find that, as well as the SCBA scattering lifetime de- 
scribed previously, a whole series of lifetimes associated 
with successive harmonics appears in the free energy. We 
define the mth lifetime r'-™' by 
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2 / ^^;ip,P>iA,,cos(m^), 



where (p is the angle between p and pi, so that t'"' co- 
incides with r and l/rtr = ^/t ~ . The free energy 
becomes F = Fq + F„ + Fj_ + Punit with 



Pn ~ TTV 



y (2^ 2. str 

•' ^ ' m=l 



PrnU.Prn, — ]^ 
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A {Pm]iP-m-l,-kk* + Pm+l,kP-m,-kk) 



F, = nvvF I d^r Str [^^^rgA (P„,P_,„_i 

m=l 

+ Pm+lP-m) — i^yA{PmP-m-l — Pm+lP-m)] , 



P± = 



F, 



-ITTVVp 



d^rStr [$^r3A(Pi -FP_i) 



+i$^A(Pi~P_i)], 
d^r Str [Po(VQ)^ -f 2iujKQ\ , (9) 



where v is the density of states, k = kx + ikyTs, $ (r) = 
U (r) VU (r), and $" ($"'") are the components of $ that 
commute (anti-commute) with A. Also Q = U {r) AU (r) 



and Do = 
rescaled P^ 



^r/2. As an intermediate step we have 



Prr 



Am) 



/r. 



Having derived the free-energy functional for the fiuc- 
tuations in (5„(r), Eq. (||), the next step is to inte- 
grate over the non-zero harmonics, represented by Pn, 
to derive a functional in terms of only U (r) . We may 
check immediately that for short-range (delta-correlated) 
disorder, the conventional unitary sigma-model is recov- 
ered: in this case, r^™) oo for m ^ 0, the P^ 
fields become infinitely massive and only the terms in 
Funit remain as required. For the case of the delta- 
correlated RMF, we find instead by explicit evaluation 
that 1/t - 1/t(") = (2m - l)/rt,., and t^"''> /t ~ 1, for 
TO ^ K^^. We then average the terms in $ with re- 
spect to the functional Pq. The leading contribution, 
containing no more than two gradient operators, comes 
from the second-order cumulant of (PI, + P_l), according 
toPo + ^;,+Pi ^ -il/2){iF,, + F^f)F,. 

To perform the averaging we apply what amounts to 
a straightforward generalization of the standard contrac- 
tion rules 1^ to include the angular harmonic space (see 
e.g. 0). The contribution of the mixed term F±F„ van- 
ishes. The terms in F± lead to a dressing of the bare 
diffusion coefficient, Dq, appearing in Punit : Dq —^D — 
VpTtr/2. Taking into account only the contributions from 
Funit and P_L corresponds to the calculation of Ref. Q| 
and gives the first line of Eq. (^). The terms in F„ in- 
duce the new term in the free energy F[U] (second line 
in Eq. (|8|)), where the dimensionless coefficient (3 can be 
written as 



f3 = 2n 



d^s 
\<so (27r)^ 



dsl 



dsl 



Tr In L(s) (10) 



Here s = k^tr/2 is a rescaled momentum variable which 
must be cut off at some upper limit set by Iq as discussed 
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above: we have sq = Ztr/(2Zo) ^ 1- The semi-infinite, 
tridiagonal matrix L in Eq. (|o|) has the entries 



{2m - l)Sm^n + is Sm+l,n + «S* S; 



m,n-t-l 



for 1 < TOjTi ^ and s = Sx + isy. We see that the 
integral in Eq. is determined only by a contribution 
of the boundary |s| = so (det L depends on s = |s| only): 



d_ 

ds 



s— |^lndetL(s~ 



(11) 



In the case where the cutoff Iq is unimportant, such as 
for sufficiently short-range disorder, this term becomes 
identically zero and the LGI is preserved as required. In 
the more general case, we see that /? is determined by 
the asymptotes of the function det L (s) , that in turn are 
determined by the range of the disorder correlations. 

For the RMF model, we have succeeded in calculating 
the determinant of L (for k ^ 0) analytically as follows. 
If we set Ln as the truncation of L to 1 < to, 7i < and 
introduce Mn{s) = detLAr(s)/((2A^ - 1)!!), we find the 
recurrence relation 



Mat = M, 



N-l 



(2N - 1){2N - 3) 



N-2, 



with the boundary conditions AIq = Mi = 1. Next we 
introduce the function f{z,s) — J2'n=o ^-'^'^ 
derive from the recurrence relation a second-order differ- 
ential equation for f{z,s). As a result we come to the 
simple formula 



N 



lim Mm{s) = cosh(s) 



(12) 



and hence /3(so) = sotanh(so) ~ sg 3> 1, which should 
be used in Eq. (^). 

We now perform a perturbative RG analysis on the free 
energy (||) (see Ref. [|]): we separate U into slow [fj) and 
fast {Uq) components, U = UUq, and expand Uq pertur- 
batively in P, where Uq = {1 - iP)'^/'^il + iP)'^''^. As 
for the conventional unitary cr-model, the diffusion coef- 
ficient D receives a negative correction at two-loop order. 
To derive the correction it is necessary to go below two 
dimensions and analytically continue the finite result to 
2D. For the free energy (^, the diffusion coefficient gains 
a further, positive correction originating from the new 
term in /?, this time at only one-loop (0(P^)) order. At 
the same time, /? is not renormalized in this order. The 
scaling function for t ~ ^/{tivD) in 2D becomes 



dt 
dlnA 



-f{l3-l) + 0{t^), 



(13) 



where A is the running momentum cutoff for the P-fields, 
with the solution 



t{u) = io [(1 + (/3 - l)tlHl/ujTtr)\ 



-1/2 



Finally we take the cutoff Iq ^ I. If we also cut off 
the range of the disorder correlations at the system size, 
L, we have k ~ l/{ppL). This gives /3 = rtr/2T ~ 
{ln{pFL)/iT^y^'^ ^ 1. We see that the resistivity de- 
creases as the frequency lo is lowered and eventually van- 
ishes, corresponding to antilocalization. 

In deriving the above scaling function it is crucial that 
the function /3(so) for sq ^ oo survives dimensional regu- 
larization below 2D. In fact, the exponential dependence 
of detL(s) on s for large s ensures that /3 remains finite 
for any dimension d > 1. For any finite range disor- 
der correlations the function det L (s) would grow alge- 
braically. Extending Eq. |ll|) to d < 2, we see that 
f3 vanishes in this case and the conventional unitary cr- 
model is recovered. 

In summary, we have considered the problem of a 2D 
electron gas in a random, static magnetic field and found 
that the resistivity vanishes in the limit of small frequen- 
cies, a behaviour which may be classified as antilocaliza- 
tion. This result may have far-reaching consequences for 
gauge field models in the theory of strongly correlated 
systems and lead to other interesting applications. 
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